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ABSTRACT 

Motivated by the irregular and little-understood morphologies of z '--^ 2 — 3 galaxies, we use non- 
parametric coefficents to quantify the morphologies of 216 galaxies which have been spectroscopically 
confirmed to lie at redshifts z — 1.8 - 3.4 in the GOODS-N field. Using measurements of ultravio- 
let (UV) and optical spectral lines, multi-band photometric data, and stellar population models we 
statistically assess possible correlations between galaxy morphology and physical observables such as 
stellar mass, star formation rate, and the strength of galaxy-scale outflows. We find evidence that 
dustier galaxies have more nebulous UV morphologies and that larger, more luminous galaxies may 
drive stronger outflows, but otherwise conclude that UV morphology is either statistically decoupled 
from the majority of physical observables or determined by too complex a combination of physical 
processes to provide characterizations with predictive power. Given the absence of strong correlations 
between UV morphology and physical parameters such as star formation rates, we are therefore un- 
able to support the hypothesis that morphologically irregular galaxies predominantly represent major 
galaxy mergers. Comparing galaxy samples, we flnd that IR-selected BzK galaxies and radio-selected 
submillimeter galaxies (SMGs) have UV morphologies similar to the optically selected sample, while 
distant red galaxies (DRGs) are more nebulous. 

Subject headings: cosmology: observations — galaxies: fundamental parameters — galaxies: high- 
redshift — galaxies: irregular — galaxies: starburst — galaxies: structure 



1. INTRODUCTION 

In the local universe the projected distribution of lu- 
minous matter within a galaxy, i.e. the morphology, of- 
ten provides a wealth of information about that galaxy's 
kinematics, rate of star formation, and recent merger his- 
tory. In the classical picture, late-type spiral galaxies 
harbor active star formation in the gas-rich arms of a 
flattened rotating disk, while early-type elliptical galax- 
ies tend to be more massive, dispersion supported, and 
quiescent systems. At high redshifts from z 2 — 3 how- 
ever the morphologies of typical galaxies are highly irreg- 
ular (Abraham et al. 1996, Kajisawa & Yamada 2001, 
Conselice et al. 2005), frequently composed of multiple 
spatially separated components, and appear to bear lit- 
tle similarity to the local Hubble-type population. It is 
uncertain whether these irregular morphologies are due 
to patchy star formation, prevalent merger activity, or 
some other physical process and consequently unknown 
whether these morphologies can (analogously to local 
galaxies) tell us anything about the star formation rate, 
mass, or stellar kinematics of galaxies at high redshifts. 

Since morphological studies are often performed at 
optical wavelengths which probe rest-frame ultraviolet 
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(UV) radiation for galaxies at redshifts z > 1, one might 
expect that the morphologies of such galaxies should ap- 
pear irregular since radiation at such wavelengths pre- 
dominantly traces emission from the brightest active star 
forming regions rather than the redder bulk of the stellar 
population (Dickinson 2000). UV emission tends to be 
patchy and irregular even for local Hubble-type galaxies 
(e.g. Gordon et al. 2004), as in the case of the local 
galaxy merger VV 114 (whose broad rest-UV absorp- 
tion line spectra suggest that it may be a local analog to 
z ~ 2 — 3 Lyman Break Galaxies; Grimes et al. 2006) 
whose near-infrared (NIR) morphology clearly shows a 
pair of interacting late-type galaxies while the rest-UV 
morphology shows only scattered clumps of emission 
(Goldader et al. 2002). However, high-redshift galaxies 
have irregular morphologies not only in the rest-UV, but 
often at rest-optical wavelengths as well (Dickinson 2000; 
Papovich et al. 2005), indicating that (in contrast to lo- 
cal galaxies) both wavelength regimes are dominated by 
emission from young starbursting components and there- 
fore that there may be some fundamental difference be- 
tween the two samples. 

One popular explanation for these multi-component, 
irregular morphologies is that they represent major 
merger systems, and that their prevalence indicates that 
the rate of major mergers was much greater at high red- 
shifts than in the local universe (e.g. Conselice et al. 
2003). Such a conclusion fits well within the framework 
of cold dark matter (CDM) theory, and may additionally 
be supported by stellar population analyses (e.g. Dickin- 
son et al. 2003) which suggest that many galaxies in the 
local universe accumulated a large fraction of their stellar 
mass at z ~ 2 — 3 as might be expected if star formation 
peaked in this epoch as a result of tidally induced collapse 



2 



Rest-UV Galaxy Morphology 



spurred by major mergers. However, the interpretation 
of a multi-component or otherwise irregular morphology 
is not always clear. In the case of VV 114, near-IR imag- 
ing (Goldader et al. 2002) indicates that all of the clumps 
of UV emission are associated with only one galaxy of 
the merger pair and that the multi-component UV mor- 
phology therefore directly reflects clumpy star formation 
rather than tracing tidally distorted features from each 
of the two galaxies. 

Building on a body of literature characterizing the 
morphologies of galaxies at redshifts z ~ 2 — 3 (e.g. 
Abraham et al. 2003; Conselice et al. 2003; Lotz et al. 
2004, 2006; Ravindranath et al. 2006) it is worthwhile 
to ask whether rest-UV morphologies correspond to 
any other physical observables such as UV/optical 
spectral line strengths (e.g. Shapley et al. 2003; Erb 
et al. 2006a), stellar population models (e.g. Shapley 
et al. 2005; Erb et al. 2006a; Reddy et al. 2006a), 
or rest-optical to IR properties (e.g. Reddy et al. 
2006a, 2006b). In this work, we use non-parametric 
coefficients to characterize the morphologies of 216 
spectroscopically confirmed galaxies in the redshift 
range z = 1.8 — 3.4, assess the statistical significance of 
correlations with spectrophotometric observables, and 
discuss the resulting physical interpretation of galaxy 
morphology. In §2 we describe our galaxy sample and 
give a basic description of the sample population. In §3 
we outline our morphological parameters, comparing our 
results to the recent studies of Conselice et al. (2003) 
and Lotz et al. (2004, 2006). Rest-frame UV spectra are 
introduced in §4, correlations between morphology and 
spectral line strength and kinematics are discussed in 
§5. In §6 and §7 we compare morphologies with stellar 
population models derived from UV to mid-IR photo- 
metric data, as well as discussing differences between 
different samples of high-redshift galaxies and AGN 
selected on the basis of various photometric criteria. 
Finally, we discuss the implications of our results for 
the physical interpretation of galaxy morphologies in §8. 
Our morphological statistics and ancillary data are made 
publically available in an electronic database located at 
|http: / / www.astro.caltech.edu/~drlaw/GOODS / 

We assume a standard ACDM cosmology in which 
Ha = 71 km s^^ Mpc"!, f],n = 0.27, and Qa = 0.73. 

2. SAMPLE SELECTION 

Our sample is drawn from rest-UV color-selected cata- 
logs of z ~ 2 — 3 star-forming galaxy candidates (Steidel 
et al. 2003, 2004; Adelberger et al. 2004) in the GOODS- 
N field. These catalogs are based on deep ground-based 
imaging, and therefore select galaxies independent of 
morphology or surface brightness since even the largest 
galaxies are nearly unresolved in these seeing-limited im- 
ages. We consider only those galaxies which have been 
spectroscopically confirmed to lie in the redshift intervals 
z = 1.8 — 2.6 or z = 2.6 — 3.4 (i.e. the peak redshift ranges 
defined by the selection functions of the color selection 
criteria, see Adelberger et al. 2004) and which exhibit no 
obvious spectroscopic signatures of active galactic nuclei. 

The redshift distribution of galaxies in our sample is 
shown in Figure ^ the z ~ 2 sample contains 150 galax- 
ies in the range 1.8 < z < 2.6 with mean z = 2.17±0.21, 
while the z ~ 3 sample contains 66 galaxies in the range 




1.5 2 2.5 3 3.5 



Fig. 1. — Distribution of galaxies with spectroscopic redshift z 
in the z ^ 2 and 2; ~ 3 galaxy samples. 



2.6 < z < 3.4 with z = 3.02 ± 0.19. 

3. MORPHOLOGIES 

Morphological parameters were determined from deep 
HST-ACS imaging obtained as part of the GOODS-N 
survey (Giavahsco et al. 2004) in F435W (S), F606W 
(y), F775W (/), and F850LP (z) bandpasses with driz- 
zled pixel scale of 50 mas pixel"^ and lOcr limiting point 
source sensitivities of 27.8, 27.8, 27.1, and 26.6 mag (AB) 
respectively (Giavalisco et al. 2004). At redshifts z ~ 2 
and z ~ 3 these bandpasses collectively probe rest-frame 
UV emission in the wavelength intervals 1500 — 3000 A 
and ~ 1000 - 2000 A respectively. 

The observed morphology is qualitatively similar 
throughout this range of wavelengths (see §3.6), and we 
therefore improve our signal-to-noise ratios per pixel by 
creating a single rest-frame UV image from a weighted 
sum of the four individual bandpasses. Weights for this 
sum are determined proportionally to the inverse vari- 
ance of the overall sky noise relative to the average num- 
ber of counts from the z ~ 2 — 3 galaxies. The UV 
composite morphologies of our 216 galaxies are shown 
in Figure |21 in order of increasing redshift and demon- 
strate a variety of morphological types ranging from sin- 
gle nucleated^ sources to extremely asymmetric sources 
with multiple nucleations and/or nebulous components. 
The "typical" galaxy has a morphology comprising one 
or more spatially distinct clumps with some degree of dif- 
fuse nebulosity, reminiscent of the IIST-5'T/S' UV mor- 
phology of the local interacting galaxy VV 114 (Goldader 
et al. 2002) which is dominated by a patchy distribu- 
tion of star formation regions. Our initial morphological 
classification groups galaxies by visual inspection on the 
basis of the apparent nucleation of their light profiles and 
the presence and number of multiple nucleated emission 

^ We adopt the term "nucleation" to qualitatively describe a 
concentrated region of flux which might naively be described as 
the "nucleus" of a given galaxy. This contrasts with the term 
"nebulosity" which we use to describe diffuse flux which is spread 
fairly uniformly over a number of pixels. 
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components. Galaxies fall within five general classes: 

1. Single strongly nucleated sources (11 sources at z ~ 
2, 9 sources at z ^ 3). 

2. Multiple strongly nucleated sources (6 sources at 
z ~ 2, 2 sources at z ~ 3). 

3. Single nucleated source accompanied by nebulosity 
(61 sources at z ~ 2, 27 sources at z ^ 3). 

4. Multiple nucleated sources accompanied by nebu- 
losity (35 sources at z ~ 2, 12 sources at z ^ 3). 

5. Nebulous emission with no strong nucleation (37 
sources at z ~ 2, 16 sources at z ~ 3). 

We seek a set of numerical parameters which will al- 
low us to effectively reproduce these intuitive divisions, 
while providing a more rigorous mathematical basis for 
the classification. The "CAS"' system of parameters has 
recently been a popular choice, characterizing galaxies 
on the basis of their concentration (C; Kent 1985, Ber- 
shady et al. 2000), asymmetry {A, Schade et al. 1995), 
and dumpiness (5, Conselice et al. 2003). However, the 
first two of these quantities are explicitly defined with 
regard to circular or elliptical apertures measured about 
a central point, which is only well defined for galaxies 
with morphologies similar to traditional elliptical or spi- 
ral galaxies, while the third quantity relies upon suitable 
choice of a smoothing scale on which dumpiness is de- 
fined. In the case of the z 2 — 3 galaxy sample, mor- 
phologies are generally so irregular (see Fig. O that they 
do not have a well-defined "center" , and the measured 
values of the CAS parameters can depend strongly upon 
the particular choice of center ^ and smoothing length 
(see Lotz, Primack, & Madau 2004, hereafter LPM04, 
for a detailed discussion). 

We therefore favor a non-parametric approach to clas- 
sification similar to that discussed by Abraham et al. 
(2003), LPM04, and Lotz et al. (2006) who define 
the gini coefficient G as a measure of the uniformity 
of the flux distribution within a source. In the follow- 
ing sections, we describe this and three additional non- 
parametric coefficients which we find effectively charac- 
terize the irregular morphologies of these z ^ 2 — 3 galax- 
ies. We note that although we considered a host of ad- 
ditional parameters in our analyses (including the Pet- 
rosian radius and a non-parametric "Petrosian area" ) , we 
found that they provided no additional information and 
therefore omit them from further discussion. 

3.1. Pixel Selection 

It is of critical importance when measuring the mor- 
phologies of faint and highly irregular galaxies to apply 
uniform selection criteria by which to assign pixels to a 
galaxy as opposed to the surrounding sky (i.e. defin- 
ing the "segmentation map" of the source) . A variety of 
criteria have been adopted in previous studies, ranging 

The asymmetry parameter A is strictly found by numerically 
searching through the image for the choice of center which min- 
imizes the value of A (see discussion by Conselice, Bershady, & 
Jangren 2000). While this relaxation technique partially mitigates 
bias arising from poor centering, it does not address the underlying 
bias present in the assumption of circular symmetry for galaxies as 
irregular as those depicted in Figure 1^ 



from complex methods based upon curve-of-growth anal- 
ysis (e.g. LPM04) to basic surface brightness selection 
(e.g. Abraham et al. 2003). 

The first of these methods, while robust to cosmolog- 
ical surface brightness dimming, can be non-trivial to 
implement in a manner consistent with non-parametric 
analysis. As outlined by LPM04, the curve-of-growth 
method calculates the elliptical Petrosian radius of a 
source (i.e. the radius from the center of the source at 
which the average fiux falls to a fixed fraction of the total 
inscribed flux; Petrosian 1976), and assigns to the seg- 
mentation map all pixels within a suitably large radius 
of this center whose fluxes are greater than the value at 
the Petrosian radius. Unfortunately, such a segmenta- 
tion map enforces elliptical Petrosian radii about a par- 
ticular center and introduces biases similar to that of the 
CAS system into the resulting morphological coefficients. 
While such a segmentation map is useful for sources with 
approximately elliptical isophotes, we find that it tends 
to fail for galaxies with multiple components or extremely 
irregular shapes since pixels at a particular "Petrosian" 
radius from an artificial center tend to include a large 
number of sky pixels, decreasing the threshhold for sur- 
face brightness selection and resulting in some fraction 
of sky pixels being allocated to the galaxy in the final 
segmentation map. While the most noticeable cases may 
be fixed by hand, this nonetheless introduces a bias as 
a function of morphological irregularity. We explore the 
effect of this bias on the measured gini coefficient in fur- 
ther detail in §3.3. 

In contrast, the second of these methods (basic surface 
brightness selection) takes no account of surface bright- 
ness dimming but is more amenable to non-parametric 
analysis. However, with the aid of our confirmed spectro- 
scopic redshifts for each galaxy in the optically-selected 
sample, we adapt this morphology-independent surface 
brightness selection technique to utilize a variable thresh- 
hold tuned to select pixels in an identical range of intrin- 
sic surface brightnesses at each redshift. 

Our segmentation map is calculated as follows: For 
each galaxy, we use our initial estimates of the position 
(based on seeing-limited UnCTZ imaging) to calculate the 
first order moment of the HST- ACS UV flux distribution 
within a 1.5 arcsecond (30 pixel) radius. A revised value 
for the center is calculated using this first order moment, 
and all pixels within a 1.5 arcsecond radius (i.e. slightly 
larger than the size of the largest galaxy in our sample, 
so the exact position of the "center" is unimportant) of 
this new center are considered as possible candidates for 
assignment to the segmentation map. While the ACS 
data product images have already been sky subtracted, 
we find that this subtraction is sometimes imperfect and 
therefore subtract off residual sky flux measured in an 
annulus of radius 1.5-2 arcseconds around the revised 
center (using a 3ct rejection algorithm to eliminate pos- 
sible contaminating flux from sources which arc nearby 
in projection). Generally, these residual sky fluxes were 
small compared to the calculated object flux. 

Once this pre-processing is complete, we assign to the 
segmentation map all pixels whose flux is at least na, 
where a is the standard deviation of pixel values in the 
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Fig. 2. — HST-ACS rest-UV morphologies of the 2^2 — 3 galaxy sample, sorted in order of increasing redshift. All panels include the 
galaxy name, redshift z, and morphological parameters size (/), gini (G), multiplicity ('f) and color dispersion (^). Images are 3 arcseconds 
on a side, oriented with North up and East to the left. Values of -2.00 for a particular morphological parameter indicate that a galaxy had 
too few pixels of suitably high surface brightness to define the parameter. Greyscale is logarithmic and chosen to emphasize the visibility 
of fainter nebulous regions, the details of high surface brightness features are thus suppressed. The 1'.'2 width of the LRIS slit is indicated 
by a dashed box in the top left panel. 



sky annulus, and where n varies with redshift as 

-3 

/ I ^ r \ 

rt = 3 
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This variable surface brightness selection compensates 
for the effects of cosmological dimming throughout our 
range of sample redshifts since (for a fixed observed band- 
pass) surface brightness scales as (l-|-z)~^ and our selec- 
tion criteria therefore include pixels of the same intrinsic 
surface brightness at all redshifts (we make the assump- 
tion that there is no intrinsic evolution of surface bright- 
ness with redshift, although see discussion in §3.2). We 
set ^max = 3.4 (i.e. the upper end of our redshift distri- 
bution) , so the value of the selection threshhold varies by 



a factor of about four across our redshift interval from 
3(7 at redshift z = 3.4 to ~ 12a at z = 1.8. We ne- 
glect the change in angular size with redshift since the 
angular diameter distance changes by only ~ 13 % over 
z ^ 1.8 — 3.4 for our assumed cosmology. The physical 
interpretation of our adopted segmentation map is illus- 
trated in Figure 13 which shows the pixels selected for a 
typical source according to 3, 5, and lOcr criteria. 

We note that by considering only flux within a 1.5 
arcsecond radius we effectively restrict our sensitivity to 
morphology on distance scales < 13 kpc at the redshift of 
our sample, similar to previous analyses (Conselice et al. 
2003a, LPM04, Lotz et al. 2006) whose selection radii 
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Figure |^ (continued) 
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range up to about 10 kpc. Based on visual analysis of 
ACS and ground-based color maps, this appears to be 
the optimal choice of distance scale to include the ma- 
jority of likely components for a particular galaxy whilst 
almost entirely excluding probable contaminants which 



appear nearby in projection. While it is likely that some 
gravitationally interacting systems extend considerably 
beyond 13 kpc, we neglect these distant components in 
our characterization of the system morphology since 1) 
Our photometric and spectroscopic data is sensitive only 
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Figure |^ (continued) 



to light at radii < 1.5 arcseconds, and 2) Close interac- 
tions are more likely than distant to produce observable 
changes in the physical state of a galaxy. 

3.2. The Size Parameter: I 



The simplest morphological parameter to define is the 
projected physical size of the source (/) seen above our 
surface brightness threshhold (we use / instead of the 
more intuitive S to avoid confusion with the "dumpi- 
ness" parameter of Conselice et al. 2003). Since the ra- 
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Figure |^ (continued) 



dius of a galaxy is not meaningful for multi-component 
systems, we use our spectroscopic redshifts to define / as 



the total projected galaxy area in square physical kpc^. 

* Given the small change in angular diameter distance across 
the redshift range of the sample there is little practical difference 
between using physical and angular sizes in our analyses. 
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Figure (continued) 




Fig. 3.— Stacked ACS rest-UV image of BX 1035 (2 = 2.238), 
overplotted with outlines of the 3, 5, and lOcr segmentation maps 
(red, green, and cyan outUnes respectively). The greyscale is linear 
in flux with the whitepoint set to lOo", the field of view is 3 X 3 
arcseconds (corresponding to a physical region 25 X 25 kpc at the 
redshift of the source). 



That is. 



/ = iV(0.05 



arcsec.2,„ , ^, ster 



pixel 



■Y[2A X 10" 



arcsec 



(2) 



where N is the total number of 50 mas x 50 mas pixels 
in the segmentation map, and Da is the angular diam- 
eter distance (in kiloparsecs) to redshift z for the as- 
sumed cosmology. This parameter makes no attempt to 
discriminate between sources on the basis of the total 
amount or relative distribution of flux within the seg- 
mentation map, and may therefore classify similarly both 
small strongly nucleated sources (which have a low num- 
ber of high flux pixels) and large yet extremely faint and 
nebulous sources (which have a low number of low flux 
pixels barely satisfying the surface brightness selection 
criteria). 

In Figure 0] we plot the distribution of / for the z ~ 
2 and z ~ 3 samples (upper left panel, black and red 
histograms respectively) . / typically ranges in value from 
less than 5 kpc^ for faint and nebulous sources to around 
10 kpc^ for the most strongly nucleated single sources, 
and up to as much as 48 kpc^ for the brightest extended 
sources. In the z ~ 2 and 2; ^ 3 samples the mean 
projected size is J = 15.0 ± 0.7 kpc^ and / = 17.4 ± 0.9 
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Fig. 4. — Relative distribution of the morphological parameters 
/ (size), G (gini), 'I' (multiplicity), and ^ (color dispersion) for the 
redshift z ~ 2 and 2 ~ 3 samples (black and red histrograms re- 
spectively). Counts are normalized by the total number of galaxies 
in each sample. 

kpc^ respectively^. This small (~ 3cr) apparent increase 
in average size from 2: ~ 2 to ^ 3 arises because the 
z ^ 3 sample does not extend as far down the luminosity 
function as the 2; ~ 2 and small faint sources are therefore 
underrepresented in the z 3 sample. If we compare 
only the fraction of the z 2 sample that overlaps the 
z ~ 3 in rest-frame UV luminosity (iieoo > 2 x lO^^L©; 
Reddy et al. 2006a) the average value of / at z ^ 2 
rises to / = 16. 7± 1.0 kpc^, which is within 1 cr of the 
z ^ 3 result. Indeed, we note that in a recent analysis 
designed to more precisely measure the size evolution of 
galaxies (as determined from their SExtractor half-light 
radii) Ferguson et al. (2004) found that the angular size 
of single-component galaxies decreases with increasing 
redshift above z '--^ 1. 

3.3. The Gini Parameter: G 

The gini coefficient (G, originally attributed to Cor- 
rado Gini [1912] and first introduced into the astro- 
nomical literature by Abraham et al. 2003) measures 
the cumulative distribution function of a population and 
may be used to distinguish between the cases for which 
flux is strongly nucleated versus uniformly nebulous (see 
LPM04 for a detailed introduction). Mathematically, G 
may be calculated as (Glasser 1962): 

1 ^ 

where X is the average flux and the Xi pixel fluxes are 
sorted in increasing order before the summation over 
all N pixels in the segmentation map. The normaliza- 
tion prefactor multiplying the summation ensures that 

^ Uncertainties quoted henceforth are uncertainties in the mean 
unless otherwise specified. 

62% of the z ~ 2 sample have well determined measurements 
of L1600, compared with 41% of the 2 ~ 3 sample. 84% of the 
2 ~ 2 sample with measured Lieoo have Lieoo > 2 X IO^^Lq. 



G takes values from zero to one inclusive, where low val- 
ues indicate a uniform flux distribution and high values 
a concentration of flux in a few pixels. 

Practically, we find that G differentiates clearly be- 
tween galaxy morphologies based on the degree of nucle- 
ation of their UV emission- sources which appear very 
nebulous generally have G < 0.2, while sources with 
strongly nucleated emission have G > 0.5 (see Fig. 
The wide range of galaxies with some combination of 
nucleated and nebulous emission components fall in a 
continuous distribution between these values. We note 
that this range of values (0 < G < 0.6) is substantially 
different from that calculated for z 2 — 4 HDF-N Ly- 
man Break galaxies (LBGs) by Lotz et al. (2004, 2006), 
who found 0.4 < G < 0.7. While some discrepancy is to 
be expected with the LPM04 results since these authors 
used rest-optical data from HST-NICMOS (although this 
difference should be minimal in light of the similarity be- 
tween rest-UV and rest-optical morphologies), we should 
expect greater similarity to the findings of Lotz et al. 
(2006) who used HST-ACS rest-UV data similar to our 
own. This difference in results appears to arise primarily 
from the adopted pixel selection method; we find that 
elliptical Petrosian selection methods similar to that of 
LPM04 tend to include more sky pixels at the "Pet- 
rosian" radius for highly irregular and nebulous objects 
(i.e. low G in our sample), which lowers the threshhold 
for assignment to the galaxy and includes more sky pixels 
in the resulting segmentation map. Including extra sky 
pixels in the map makes the genuine galaxy pixels appear 
to be comparatively more rich in flux, artificially increas- 
ing the value of G calculated for the sources. Applying a 
segmentation map similar to that of LPM04, we find that 
the gini coefficient calculated for our highest G sources 
remains relatively unchanged, while our low gini sources 
increase their values of G considerably, artificially com- 
pressing the range of values to 0.4 < G < 0.6, closely 
mimicking the range of values presented by LPM04. We 
conclude that while the LPM04 values are accurate in 
the sense that a G = 0.4 (in their calculation) galaxy is 
more nebulous than a G = 0.6, the noise introduced by 
this compression in dynamic range severely hampers the 
discriminating power of the gini coefficient. 

A histogram of values of G is plotted in Figure 01 and 
suggests that the mean value of the gini coefficient ap- 
pears to change from G = 0.27 ± 0.02 at z ~ 2 to 
G — 0.35 ± 0.01 at z ~ 3. To some extent this may 
be due to the underrepresentation of faint, nebulous ob- 
jects in the 2 ~ 3 sample (see §3.2), but even restricting 
the sample to the same range in intrinsic UV luminos- 
ity the difference between the two samples is of order 3cr 
suggesting that the UV emission from objects at z 3 
is genuinely slightly more nucleated than at z ~ 2, con- 
sistant with the finding of Ferguson et al. (2004) who 
demonstrated that apparent galaxy size at constant lu- 
minosity decreases with increasing redshift. 

We note that to first order, G does not distinguish 
between sources based upon their number of nucleated 
components- i.e. a galaxy with two or more apparent 
nucleations has a value of G nearly identical to a galaxy 
with only one, so long as the cumulative distribution of 
light is similar. That is, the exact spatial distribution of 
flux is irrelevant to G, which is instead sensitive to the 
overall curve of growth of the total flux. 
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3.4. The Multiplicity Parameter: ^' 

Our third classification parameter (5") is designed to 
discriminate between sources based on how many ap- 
parent components the light distribution is broken into, 
i.e. how "multiple" the source appears. This parameter 
is similar to both the asymmetry parameter A (Schade 
et al. 1995) and the second order moment of the 20% 
brightest pixels (M20, LPM04) in that it is sensitive to 
the presence and distribution of multiple clumps of flux. 
However, both A and M20 have their limitations: A de- 
pends strongly on the assumption of overall circular sym- 
metry about some central point for each galaxy, while 
M20 is normalized by the moment of the segmentation 
map to remove the effect of overall galaxy size, which 
unfortunately results in a limited dynamic range since 
the segmentation map and the 20% brightest pixels of- 
ten have a similar spatial distribution. In contrast, 5* 
is defined in a manner that requires neither a center of 
symmetry nor a conventional normalization. 

Using the observed flux distribution as a proxy for 
"mass" , we calculate the "potential energy" of the light 
distribution projected into our line of sight as 
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where the summation of the pixel fluxes Xi and Xj runs 
over all TV pixels in the segmentation map, and is the 
distance (in pixels) between pixels i and j. This value 
is normalized by that which would be achieved with the 
most compact possible re-arrangement of the flux pix- 
els, i.e. the configuration which would require the most 
"work" to pull apart. We re-arrange the physical po- 
sitions of the N pixels of the segmentation map in a 
circular configuration with the brightest pixel in the cen- 
ter and with pixel flux decreasing outwards with radius. 
Calling the distance between pixels i and j in this re- 
arranged map r'^j, the projected potential energy of this 
compact light distribution is 



N 



N 



V'compact — ^ ^ ^ ^ 

i=l j = l.j^i 



XiXi 



(5) 



The multiplicity coefficient is then defined logarithmi- 
cally as the ratio of these two quantities: 



* 100 X logi 
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(6) 



factual 

As shown in Figure^l ^E" ranges from up to about 30. 
Typically, = — 2 for single well-nucleated sources, 
\1/ ~ 5 for sources which are beginning to show a sec- 
ond component in addition to the main nucleation, and 

~ 10 for strong double-nucleation systems. At values 
of 5* > 10, ^' increases as the number and separation 
distance of nucleations increase, but numerical distinc- 
tions become less obvious until at 5' > 25 the major- 
ity of sources are low G nebulous sources for which 5" 
breaks down as a useful statistic (see §3.6). There is 
therefore some degree of correlation between and the 
gini coefficient G (Fig. I^J since sources with very neb- 
ulous emission (i.e. small G) tend to have very spread 
out fiux distributions (i.e. large ^). However, this is 
of secondary importance to the largely orthogonal clas- 
sification permitted by these two parameters; while G 
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Fig. 5. — Distribution of G (gini) versus ^ (multiplicity) for 
galaxies in the z 2 (black squares) and 2; 3 (red triangles) 
samples. Values are correlated in the sense that faint, nebulous 
galaxies tend to be broken into a greater number of non-contiguous 
components. 



is most sensitive to the overall curve of growth of the 
flux distribution and fairly insensitive to the number of 
nucleated components, \1/ is sensitive to the number of 
bright components and comparatively insensitive to the 
overall curve of growth. We note that since ^E" is a flux- 
weighted statistic it is more sensitive to bright pixels, 
and therefore a bright central source with an extremely 
faint secondary source will tend to have low while the 
same central source with a bright secondary source will 
tend to have higher . 

3.5. The Color Dispersion: ^ 

Finally, we capitalize on the available multi- wavelength 
HST-ACS data by defining the color dispersion parame- 
ter ^ (Papovich et al. 2003) , which quantifies the degree 
of morphological difference between two bandpasses. Ap- 
plied between rest-optical and UV wavelengths, f can 
measure the difference in spatial distribution between 
stellar populations of differing ages, convolved with vari- 
ations in the dust distribution and resulting extinction 
(Papovich et al. 2003). Using rest-UV data alone ^ is a 
less powerful statistic, but nevertheless potentially infor- 
mative. 

We calculate as 



J2{h~aIv-f3)'-j:iB. 



aBvf 



(7) 



where Iv,Iz are the object pixel fiuxes in V and z band- 
passes respectively^^, By.Bz the background sky flux 
in the bandpasses, a is the flux ratio between the two 
bands, and (3 represents the different in background lev- 
els between the two bands. In brief, the first term in 
the numerator represents the summed square difference 

While use of the B and z bands would provide a greater 
wavlength baseline for morphological differences, we use V instead 
since B is blanketed by absorption from the Lyman a forest for 
galaxies in the upper end of our redshift range. 
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Fig. 6.— HST-ACS images of BX 1157 and BM 1139 in F606W 
(Panels a and b respectively) and in F850LP (Panels c and d). Both 
galaxies have comparable redshifts (z = 2.08 and 1.92), gini (G = 
0.291, 0.321), and multiplicity (I- = 3.23, 4.94) in the summed 
rest-UV image. Both galaxies show similar decrease in signal-to- 
noise ratio from F606W to F850LP, but have widely variant color 
dispersions (^ = 0.023,0.112). Image size and orientation are the 
same as in Figure 1^ 



in pixel fluxes between the bands, the first term in the 
denominator the summed square total object flux (for 
normalization), and the second terms in both numera- 
tor and denominator represent corrections to statistically 
eliminate contributions to ^ from the natural background 
sky variance. Further details regarding the definition of 
this statistic are given by Papovich et al. (2003). 

As shown in Figure^ typical values range from ^ ^ 0.0 
to 0.15, with f = 0.057 ± 0.006 at z ~ 2 and f = 
0.061 ± 0.011 at z 3 (i.e. consistent with a constant 
value to well within the uncertainty). We note that this 
range of values is larger than that of ^ 0.0 — 0.05 found 
by Papovich et al. (2005) for a sample of z ^ 2.3 galax- 
ies measured between rest-UV and optical bandpasses. 
The origin of this discrepancy is uncertain, although it is 
likely that different galaxy samples, segmentation maps, 
and bandpasses all contribute. By and large, it appears 
that our calculated ^ is dominated mainly by scatter 
rather than by genuine differences between the appar- 
ent morphology in V and z bands. Figure IHl shows that 
for two galaxies with ^ — 0.023 and ^ = 0.112 the most 
obvious trend from V to z bands is an overall decrease 
in the signal-to-noise ratio. Given the similarity of the 
effect in both cases, it is not readily apparent why the 
two galaxies should have such different color dispersions. 

Curiously, while ^ has no obvious correlation with a 
visible difference in morphology between bands, it does 
correlate with the Gini coefficient G of a source (as shown 
in Fig. Cj) in the sense that ^ is (on average) slightly 
larger in the most nucleated sources. This is likely a 
consequence of the flux-weighting of ^- high flux pixels 
have greater absolute variation between bandpasses and 
tend to dominate the sum in Equation 7, increasing ^ for 
high G sources where f is dominated by variation from 
a few bright pixels rather than averaging over a larger 



0.4 - 



0.2 



0,2 



0,4 



0,6 



Fig. 7. — Distribution of G (gini) versus g (color dispersion) 
for galaxies in the 2; ~ 2 (black squares) and z ~ 3 (red trian- 
gles) samples. Note that a few extreme outliers from this trend lie 
outside the visible region of this plot to better illustrate the main 
correlation. 



number of lower- variation pixels for low G sources. In 
conclusion, we caution that while ^ is ideally a useful 
statistic, in the present case it may be too erratic to 
provide a great deal of information. 

3.6. Robustness of the Parameters 

We test the robustness of these parameters to cosmo- 
logical distance by selecting one representative galaxy for 
each of the five general morphological categories (defined 
previously in this section), artificially redshifting them 
through the range z — 1.8 — 3.4, and measuring the re- 
sulting morphologies using our variable threshhold pixel 
selection technique. As illustrated by Figure |S1 in most 
cases we would measure consistent values for the morpho- 
logical parameters for a given galaxy if it were located 
at any redshift throughout our sample. This constancy 
fails however for the most faint and nebulous of objects 
(G < 0.15, green dot-dashed line in Figure |H1), for which 
the multiplicity ^E" and the color dispersion ^ can vary 
subtantially because the number of pixels in the segmen- 
tation map is small and random variations in noise can 
drastically affect both the locations of selected pixels (to 
which is particularly sensitive) and the residual color 
(to which ^ is sensitive). Such objects represent less than 
20% of the sample however, and the importance of this 
effect on ^ may be mitigated by noting that is only 
well-defined up to ~ 25. We also note a slight decrease 
in / with redshift to z = 3 , followed by a smaller upturn 
at z > 3.3. The first of these effects is consistent with 
our neglect of the change in angular size with redshift, 
and the second with the inclusion of a small number of 
sky pixels in the segmentation map when the selection 
threshhold is pushed down to 3fT. Neither effect is large 
enough to noticeably impact our analysis. 

Having demonstrated the uniformity of the selection 
technique, we consider what (if any) information is lost 
for galaxies in the lower end of the redshift range by ef- 
fectively restricting our analysis to those pixels brighter 
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Fig. 8. — Figure illustrates the robustness of morphological parameters to redshift when using the variable surface brightness method to 
define the segmentation map. Colored lines indicate the trends found by artificially redshifting one representative galaxy from each of the 
five general morphological classes through the range z = 1.8 — 3.4 before performing morphological analysis. These five classes (defined in 
§3) and their representative galaxies are: Single nucleated source (BX 1040: solid red line), multiple nucleated source (BX 1630: dotted 
magenta line), single nucleated source with nebulosity (BX 1297: short-dashed dark blue line), multiple nucleated source with nebulosity 
(BX 1035: long-dashed cyan line), and purely nebulous source (BX 1169: dot-dashed green line). 



than about lOcr above the background sky noise. In Fig- 
urelHlwe plot morphological parameters calculated for the 
same five galaxies as before for a range of surface bright- 
ness selection threshholds. Most obviously, / decreases 
rapidly with increasing brightness threshhold since corre- 
spondingly fewer pixels are included in the segmentation 
map. In addition, G declines noticeably as the lower- 
flux population is gradually omitted from the map. In 
contrast, and ^ remain relatively constant throughout 
the range of selection threshholds considered, with the 
exception of purely nebulous sources (dot-dashed green 
line) for which we have previously noted the instabil- 
ity of these two parameters. Most importantly however, 
if we neglect purely nebulous sources we note that the 
same relative information is preserved at all brightness 
threshholds- no matter what the threshhold it is equally 



possible to distinguish between the galaxy types despite 
the overall trends. That is, since the lines do not cross 
(again with the exception of the green dot-dashed line 
and ^ at very low threshholds) we may be confident that 
we are not discarding information by using a ~ lOtr seg- 
mentation map at lower redshifts. 

We now note a few caveats to this analysis which 
should be borne in mind. First, by using a fixed bandpass 
we have measured the morphologies of z ~ 2 and z ~ 3 
galaxies at slightly different rest wavelengths. However, 
this effect should be negligible since our morphological 
parameters do not change substantially using data from 
the B, V, or / bandpasses (although they can change 
somewhat in the z bandpass due to the lower signal-to- 
noise ratio for all detections at near-IR wavelengths). Pa- 
rameters calculated in each of these bands correlate very 
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Fig. 9. — As Figure |S] but plots the calculated values of morphological parameters for the five classes of galaxies as a function of the 
adopted surface brightness selection threshhold in multiples of the background sky noise cr. Note for comparison that segmentation maps 
corresponding to the 3, 5, and lOcr selection regions for BX1035 (i.e. long-dashed cyan line) are shown in Figure 1^ 



strongly with one another, differences being dominated 
by the hmiting surface brightness reached in each band. 
Given that our sample consists of actively star-forming 
galaxies, this is perhaps not surprising since young mas- 
sive stars tend to dominate the SED out to at least rest- 
frame blue wavelengths. Second, we have neglected sur- 
face brightness evolution throughout our sample. Since 
our primary goal is to investigate the physical meaning of 
morphology rather than providing a detailed comparison 
between samples at different redshifts we choose not to 
introduce uncertain corrections for such effects into our 
analysis. Third, it is possible that some regions of flux 
included in our segmentation maps are due to chance 
alignments of some object at a very different redshift 
than the target galaxy. However, we judge this to be 
an unlikely major source of contamination to our sample 
since 1) the various pieces of the galaxies have consistent 
colors, and 2) the total number density of optically se- 



lected z 2 — 3 galaxies is ~ 10 arcmin"^ to TZ = 25.5 
(Steidel et al. 2004), giving a low probability for chance 
supposition of two unrelated galaxies within 1.5 arcsec- 
onds. Finally, we have made no correction for the effects 
of galaxy inclination to the line of sight. While it is 
conceivable, and indeed likely, that some of the variance 
in apparent morphology may be due simply to different 
orientations relative to the line of sight, our knowledge 
of these galaxies is at present insufficient to allow us to 
compensate for such effects in any meaningful way. 

Of the four parameters which we have introduced, the 
G — ^I* classification scheme best reproduces the morpho- 
logical trends apparent to the eye. This system is illus- 
trated in Figure [TUl which shows a representative set of 
galaxies from our sample and demonstrates how these 
two parameters serve to distinguish galaxies based on 
their degree of nucleation and number of components. 
Since surface brightness dimming and bandshifting make 



14 



Rest-UV Galaxy Morphology 



it difficult to associate the morphologies of z ~ 2 — 3 
sources with galaxies in the local universe we avoid such 
direct, and potentially misleading comparisons. We de- 
scribe the sense of our classification parameters by not- 
ing, however, that early-type galaxies typically are larger 
(i.e. larger /), more concentrated (larger G), less multi- 
ple (lower and have lower color dispersion (^) than 
their late- type counterparts which are instead dominated 
by emission from multiple scattered knots of star forma- 
tion. A more thorough overview of the relation between 
local and high-redshift galaxies is given by LPM04. 

4. OVERVIEW OF REST-UV SPECTRA 

4.1. Spectral Processing 

As part of an ongoing Keck LRIS-B spectroscopic sur- 
vey (Steidel et al. 2004) we have compiled UV spectra of 
our 216 target galaxies in the GOODS-N field, and here 
explore how these spectra correspond to the morphology 
of their host galaxies. Unfortunately, spectra of individ- 
ual galaxies are rarely of sufficient quality to accurately 
measure the strengths of their emission/absorption com- 
ponents, and we therefore divide the spectra into five bins 
according to each of our four morphological parameters 
and measure the strengths of features in stacks of the 
spectra within these bins. We find that five bins gives a 
suitably large number of bins across which spectroscopic 
trends may be assessed while still producing reasonable 
quality stacked spectra (each comprised of 30 galaxies at 
2 ~ 2 and 13 galaxies at z ~ 3). 

Our spectroscopic combination method is similar to 
that described by Shapley et al. (2003). Before stacking, 
individual sky-subtracted spectra were ffux-calibrated 
and shifted to the systemic rest frame using the prescrip- 
tions given by Adelberger et al. (2005), then resampled 
to a common dispersion of lA pixel" ^ and rescaled to 
a common mode in the range AA1250 - 1500 A. Spec- 
tra were manually cropped to eliminate overly noisy seg- 
ments due to falloff of the LRIS-B blue-side efficiency 
and of the red-side dichroic transmission. The resulting 
spectra were averaged using a min/max rejection of hot 
pixels, cosmic rays, and bad sky subtraction events; three 
high and low values were rejected at each dispersion co- 
ordinate (corresponding to a rejection of the high/low 
10% for the z ~ 2 sample). Continuum normalization 
was performed by iteratively fitting a spline function to 
the stacked spectrum using wavelength intervals selected 
to be free of strong interstellar features, giving results 
consistent with the normalization method described by 
Rix et al. (2004). 

4.2. Key Spectral Features. 

The rest-frame ultraviolet (AAIOOO — 1800 A) spectra 
of rest-UV color selected galaxies are typified (Fig. Ill|l 
by redshifted Lyman a (Lya) emission superposed on 
a broad blueshifted resonant absorption trough, and by 
strong velocity-broadened absorption lines due to ener- 
getic galaxy-scale outflows (Shapley et al. 2003). The 
strength of Lya emission is originally governed by the 
rate of star formation and initial mass function (IMF) 
in galactic Hii regions, although this raw flux is signifi- 
cantly modified by resonant scattering from interstellar 
H I, generally culminating in its eventual absorption by 
interstellar dust and subsequent reradiation in the in- 
frared. The observed Lya emission strength is therefore 



a complex function of neutral hydrogen column density, 
dust fraction, and geometric/kinematic projection of the 
outflowing ISM on the line of sight. Typically, observed 
Lya emission is redshifted by roughly 400 km s^"'^ rela- 
tive to the systemic velocity (Shapley et al. 2003, Steidel 
et al. 2004), corresponding to scattering from the back 
side of outflowing ISM, from which the Doppler shift is 
suitable to take the photon off of resonance and permit 
it to escape the galaxy. 

Other notable features include strong, optically 
thick absorption lines from low-ionization species (e.g. 
Sin A1192, Sin A1260, Oi-I-Sn A1303, Cn A1334, 
Sin A1526, Fen A1608, and Aln A1671) arising from 
outflowing neutral ISM. As shown in Figure ^2 these 
absorption lines are blueshifted from the systemic red- 
shift, corresponding to absorption from outflowing gas 
seen from the nearby side of a galaxy-wide outflow — 
their full widths can reach 1000 km s""'^ (Pettini et 
al. 2002; Shapley et al. 2003). We note that while 
Oi-I-Sin A1303 (a blend of the low ionization species 
Oi A1302 and Sin A1304) appears to be blended with 
a third component around A1297 (resulting in a shal- 
lower blue-side slope for the composite line) it is gen- 
erally possible to deblend this additional component to 
measure Oi+Sin A1303 alone. This broad A1297A fea- 
ture is likely itself a blend of the stellar photospheric 
lines Cm A1296.33, Sim A1294.54, Sim A1296.73, and 
Sim A1298.93 (Tremonti et al. 2006), variations in the 
relative strengths of which can shift the apparent cen- 
troid of the A1297 A blend from A = 1295 to 1298 A. 

Also apparent are absorption lines due to higher- 
ionization species, including the Siiv AA1393, 1402 dou- 
blet and Civ absorption around A1549A. Although all 
of these lines are blended with stellar features, C iv is 
a particularly complex blend of interstellar absorption 
lines at A1548A and 1550 A, combined with a P-Cygni 
component from the winds of the most luminous O and B 
stars. We neglect the numerous additional features due 
to nebular, fine-structure, and stellar atmosphere transi- 
tions (e.g. Cm A1176; see Shapley et al. 2003 for a fur- 
ther list), since they are not generally detected to high 
statistical significance after the sample has been divided 
amongst five morphological bins. We refer the reader to 
Shapley et al. (2003) and references therein for further 
discussion of the physical interpretation of the rest-UV 
spectra. 

4.3. Equivalent Widths and Uncertainties 

The equivalent widths of all absorption features were 
integrated non-parametrically relative to the normal- 
ized continuum. Generally, features shortward of Lya 
and longward of A1500 A are considerably noisier than 
those between these wavelengths, and we therefore de- 
fine average "low-ionization" and "high-ionization" ab- 
sorption line strengths (M^lis and Whis respectively) as 
the weighted means of the shorter wavelength transitions 
Sin A1260, Oi-^Sin A1303, and Cn A1334 (for Wlis) 
and Si IV A1393A and A1402A (for I^ms)- 

The total equivalent width of the Lya feature (i.e. ab- 
sorption plus emission) tends to be strongly affected by 
the noisiness of the absorption trough, so we character- 
ize the strength of the emission component alone. The 
equivalent width of emission (WLya) is determined by di- 
viding the total flux in the emission component by the 
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Fig. 10. — As Figure 12] HST-ACS rest-UV morphologies of redshift z ~ 2 galaxies classified according to gini (G) and multiplicity 
(^) parameters. Horizontal and vertical bin ranges are chosen to divide the sample into quartiles. Increasing values of G correspond to 
increasing nucleation of source emission, and increasing values of correspond to increasing number of components. 



continuum level which would be present in the absence 
of any absorption trough. 

It is possible to assess the statistical significance of 
possible deviations from a constant value by using the 
statistic 

^ \ ^'sample / 

where Xi is the measured equivalent width for a given 
bin, X is the mean equivalent width among all five bins, 
and cTsampic is the square root of the natural variance 
among samples of randomly binned spectra. We find 
that asampic is considerably larger than the uncertainty 
in the mean equivalent width for a given stack (which 
belies the actual variance observed between randomly- 
drawn samples), and therefore adopt it as a conservative 
means of ensuring that any possible trends are more sig- 

Thirty samples was sufficient to converge o-gamplc to within 

5%. 



nificant than would be likely to occur randomly. As such, 
the significance of a particular value of is evaluated 
using 5 — 1 = 4 degrees of freedom to give the likeli- 
hood P that deviations of the five measurements from a 
constant average value (whether in the form of a mono- 
tonic trend or a single bin whose equivalent width varies 
greatly from the average with respect to the expected 
variance) are greater than that expected from a random 
binning of the sample. Applying this method to our mea- 
surements from randomly selected quintiles of galaxies, 
we determine that a threshhold of P > 90% suffices to 
cull apparent associations due to random variance. 

5. THE RELATION OF REST-UV MORPHOLOGIES TO 
SPECTRA 

As described in §4.1, we bin the galaxy sample sep- 
arately into quintiles according to each of the morpho- 
logical parameters /, G, and ^. Bin divisions are 
determined so that each bin contains an identical num- 
ber of galaxies (i.e. 30 each for the z ~ 2 sample, and 13 
each for the z 3 sample), precise ranges are given in 
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Fig. 11. — Stacked continuum- normalized spectra are plotted for all 150 2 ~ 2 and 66 2 ~ 3 galaxies (black/red lines respectively). 
Labels and dotted lines indicate the fiducial locations of atomic transitions corresponding to major spectroscopic features. Note that Lya 
emission is considerably stronger in the 2; ~ 3 sample than the 2 ~ 2, and peaks at a relative normalized flux of 3.95 (not shown for clarity). 



Tabled In all cases our spectra are broadly consistent 
with the diffuse Ught spectrum of local starburst galaxies 
(e.g. Chandar et al. 2005) indicating that largely simi- 
lar processes likely dominate the UV light output of all 
morphological types in our sample. Despite this general 
similarity the spectra show some variation with morphol- 
ogy as demonstrated in Figure^] which we proceed to 
discuss in detail. 

We first note however the caveat that the position an- 
gles of the slits used to obtain our UV spectra have not 
been chosen to correspond to the major axis of each of 
the galaxies, and it is therefore possible for a few of the 
most widely separated sources with multiple components 



that the UV spectra may represent only one of the com- 
ponents. Given the seeing limited nature of the spec- 
troscopic observations and the fact that the typical size 
of our targets 1 arcsecond) is less than the width of 
the LRIS-B slit (1'.'2) we doubt that this has a consider- 
able effect. More likely perhaps is the probability that 
spectra are dominated by light from the bright nucleated 
regions of our sources and may not be expected to show 
any difference between nucleated and nucleated -I- nebu- 
lous sources if the spectrum of the nucleated regions are 
similar in each case. 

5.1. Interstellar Absorption Lines 
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Fig. 12. — LycK emission {Wi^ya) equivalent width, low- and high-ionization interstellar absorption equivalent width (Wlis and H^his 
respectively), and the kinematic offsets between emission and absorption lines (Atiem-abs) are shown measured from quintiles in the 
morphological parameters I (size), G (nucleation), 't (multiplicity), and ^ (color dispersion). Black and red points represent data from the 
z ^ 2 and z ~ 3 samples respectively, shown with error bars representing the standard deviation among measurements made from randomly 
binned samples. Red and black dotted lines indicate the values measured for the complete z ~ 2 and 2 ~ 3 stacks respectively, while the 
black and red percentages given in each panel indicate the probability P that the data points shown are statistically inconsistent with 
a constant value. Generally, values of P > 90% are indicative of significant deviation- note particularly the linear trends of Lya emission 
increasing for more nucleated sources (i.e. those with higher G) and low-ionization absorption strength increasing for larger sources (i.e. 
those with greater /). 



As illustrated in Figure ^] (middle rows), we find 
that the strength of interstellar absorption lines (Wlis 
and Whis) is largely uncorrelated with UV morphology, 
although there appears to be a statistically significant 
(P = 92%) association of the strength of low-ionization 
species with galaxy size /. Indeed, the data may be con- 
sistent with a trend that larger galaxies tend to have 
stronger interstellar absorption lines. Such a trend is 
most apparent for the z ~ 2 sample; while the z ^ 3 
sample suggests a similar trend the greater statistical 
uncertainties mitigate its significance. 

We explore this possible trend in greater detail for the 
z ~ 2 sample by stacking our spectra in two dimen- 



sions to explore the variation of line strengths across a 
plane of two parameters smoothed by a variable-width 
kernel. We stack the spectrum of each galaxy with that 
of its 10 nearest neighbors in the I — G plane to dis- 
tinguish large nebulous objects from large yet nucleated 
objects, with inter-point "distances" along each axis nor- 
malized by the dynamic range of values along the axis. 
In Figure [TSI we plot / versus G for z ~ 2 galaxies, with 
point size corresponding to the strength of M^lis- On 
the whole, the distribution of galaxies with strong and 
weak low ionization lines (i.e. large and small points) is 
quite similar, except for the lower left corner represent- 
ing the faintest and most nebulous galaxies, for which 
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TABLE 1 

QUINTILE BINS FOR MORPHOLOGICAL PARAMETERS. 



Parameter 




X2 


xs 


X4 


^5 


X6 


I {z^ 


.2) 


0.181 


7.473 


12.540 


15.977 


20.108 


48.163 


I (z^ 


-3) 


4.200 


11.559 


13.952 


17.991 


21.407 


45.181 


G {Z r 


-2) 


0.034 


0.172 


0.245 


0.321 


0.382 


0.568 


G {Z r 


-3) 


0.105 


0.265 


0.321 


0.385 


0.459 


0.570 


( Z r 


-2) 


-0.199 


1.667 


3.871 


6.151 


10.677 


34.062 


"I- {Z r 


-3) 


0.249 


2.426 


5.236 


8.849 


12.599 


22.795 




.2) 


-0.115 


0.015 


0.033 


0.049 


0.076 


0.507 




-3) 


-0.095 


0.008 


0.040 


0.057 


0.078 


0.500 



ously distinct class of galaxies, or have any obvious rea- 
son for occupying the region of this plot which it does. 
We suggest the possibility that there could be a popu- 
lation of galaxies which have relatively weak outflowing 
components which are somehow distinguishable in this 
plot, but caution that random variation combined with 
the kernel smoothing technique might also be responsible 
for this enigmatic feature. 

There is one additional significant deviation from uni- 
formity for the color dispersion ^ and the low-ionization 
absorption line strength (at a significance level of > 99% 
for the z ~ 2 sample). However, inspection of Figure [T^ 
offers no clear explanation of the nature of this associa- 
tion: objects in both the smallest and largest bins of ^ 
have weaker absorption lines than objects in the interme- 
diate three bins, and there hence appears to be no par- 
ticular linear correlation between the variables. While 
this deviation may be a genuine product of physical phe- 
nomena its interpretation is unclear, particularly given 
the uncertain meaning of ^. 

5.2. Ly a Emission 

As indicated by Figure [T^ fbottom row) there are rela- 
tively many possible associations between Lya emission 
strength and galaxy morphology, the simplest of which 
to interpret is the positive correlation between Lya and 
the gini coefficient G. These two parameters show a clear 
linear trend of increasing emission strength with increas- 
ing nucleation; for the z ~ 2 sample the most strongly 
nucleated sources show roughly twice the emission seen 
in the most nebulous (see also Fig. I14|l . and for the z ^ 3 
sample an even stronger trend is possible (although less 
well defined). It is unlikely that this trend is due to 
substantial variation of the composition or velocity of in- 
terstellar gas since such variations should also affect the 
strengths of absorption lines, which remain statistically 
constant across all values of G (see Figs. 1121 and lT^ . A 
more plausible explanation may be that dust might be 
present in greater quantities in nebulous low-G sources, 
causing greater attenuation of Lya photons. Under this 
hypothesis, it may be simply the presence of more or less 
dust which determines both the observed degree of UV 
nucleation and the strength of resonant Lya emission. 
We explore this hypothesis further in §6. 

We map the Lya trend in greater detail in Figure El 
which compares the strength of the association between 
Lya and G with the previously known association be- 
tween Lya and stellar mass M^, (i.e. that galaxies with 
higher stellar masses have lower average Lya emission 
strength, Erb et al. 2006b) Figure [T31 suggests that 
the overall association between G and Lya is genuine, 
but that emission strength in fact peaks for galaxies with 
G ~ 0.4 rather than for the few galaxies with G > 0.5. 
Since the galaxies with G > 0.5 represent the extreme 
of nucleation and have much weaker Lya emission than 
expected based on the majority of the galaxy sample we 
posit that these few galaxies may be somewhow distinct 
from the rest of the sample. Alternatively, the G ~ 0.4 
galaxies may represent a particularly dust-free popula- 
tion from which it is possible to see both star forming 

We caution that despite the fact that G is correlated with 
Wi^yc, and Wi^yc, with M* , G itself is not correlated with M* (see 
§6), indicating that correlations are not necessarily commutative. 



^Successive values of Xi denote boundary divisions between galaxies 
binned into equal-size quintiles according to each of the morphological 
parameters (without regard to the other three parameters). 
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Fig. 13. — Low-ionization absorption line strengths for z ^ 2 
galaxies are plotted (plusses) according to their location in the 
size — gini (I — G) plane, spectra of each galaxy have been stacked 
with those of their ten nearest neighbors in this plane. Point sizes 
correspond linearly to the value of Wlis- Note the overabundance 
of weak-absorption sources (smaller symbols) in the lower left cor- 
ner of the plot, and in the strip running from G, I = 0.2, 30 kpc^ 
to 0.6, 10 kpc^. 



there is an overabundance of small points (i.e. sources 
with weak Wlis)- Given the similar overall distribution 
of line strengths, we conclude that the association be- 
tween / and low-ionization absorption strength is due to 
the over-representation of weaker-line sources among the 
faint and nebulous galaxy sample rather than to an over- 
all trend. Although this may indicate a genuine physical 
characteristic, we note that it is this class of low surface- 
brightness galaxies for which spectra are typically of the 
poorest quality and absorption line measurements least 
reliable. It is interesting, however, to note that the high 
ionization absorption lines do not likewise appear weaker 
in this class of galaxies (see Fig. I12|l as might be expected 
were the apparent decline in low-ionization line strength 
due to spectrum quality alone. 

Figure 1 11^1 also introduces another curiosity, namely 
the overabundance of weak-line objects in a line run- 
ning through the plot from G,I — 0.2, 30 kpc^ to G, / = 
0.6, 10 kpc^. This band does not correspond to an obvi- 
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Fig. 14. — Stacked, continuum-normalized spectra for 2 ~ 2 
galaxies representing three bins in the morphological nucleation 
parameter G. The low G nebulous galaxy bin is represented by 
the black line, the intermediate bin by the cyan line, and strongly 
nucleated sample by the magenta line. Spectra have been smoothed 
by a three pixel boxcar filter and normalized consistent with the 
prescriptions of Rix et al. (2004). Labels and dotted lines indicate 
the fiducial locations of major spectral features. 
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Fig. 15. — Lyo emission line strengths for z 2 galaxies are 
plotted (plusses) according to their location in the stellar mass — 
gini (M* — G) plane, spectra have been stacked with their ten 
nearest neighbors in this plane as described in §5.1. Point sizes 



correspond linearly to the value of Wi^y 



Note the concentration 



of sources with strongest Lyo emission around G ~ 0.4 and A/* 
1.5 X 10^°Mq. 



regions and their surrounding material. 

We also observe also that Wlyq deviates from statisti- 
cal uniformity when binned according to multiplicity 5* 
(for both z ~ 2 and z ~ 3 samples), apparent size / (for 
the 2; ^ 3 sample), and color dispersion ^ (for the z ~ 3 
sample). In the first case, there is less of an obvious 
trend in Figure [T^ than an apparently random scatter of 
points, with particularly discrepant points found in both 
the highest and lowest bins in In the latter two cases. 



W^Lya possibly declines with / and increases with ^ for the 
z ~ 3 sample, but these trends are absent from the z 2 
spectra. It is possible that these apparent deviations are 
due to the residual correlation between our morphologi- 
cal parameters (in which case the sense of a trend would 
be a complicated projection of the G — Whya correlation 
onto the density distribution in morphology space), or 
simply that these deviations are telling us of some trend 
or population of galaxies which is not well distinguished 
using our chosen morphological parameters (or indeed 
any of the many others which we have explored). How- 
ever, we are inclined (despite our conservative estimates 
of the uncertainty) to ascribe at least those deviations 
only visible in the z ~ 3 spectra to statistical variance 
since similar trends do not appear in the much higher- 
quality composite spectra from the z ~ 2 galaxy sample. 

5.3. Kinematic Offsets 

Since our stacked spectra have all been shifted to the 
systemic rest frame (i.e. that of the stars and [H 11] re- 
gions), it is possible to calculate the kinematic offset 
Awcm-abs in cach stack between Lya emission and the 
average of the low-ionization absorption lines Sill A1260, 
Oi+Siii A1303, and C II A1334. For the z ~ 2 and z ~ 3 
samples respectively we find that Aiiom-abs = 453 and 
594 km s~^, roughly similar to the Shapley et al. (2003) 
value of Awcm-abs — 510 km s^^ found for a large stack 
of 794 z 3 LBG spectra While our relatively small 
sample of 66 spectra at z 3 is too noisy to distinguish 
any significant trends with morphology, for the z ~ 2 
sample Aiicm-abs is convincingly correlated with both 
the size parameter / and the gini coefficient G (see Fig. 
1121 top panels). 

In particular, the offset velocity increases nearly lin- 
early from ~ 300 km s~^ for the smallest galaxies with 
J ~ 5 kpc^ to near 600 km s~^ for the largest with / ~ 20 
kpc^. The differences are similarly pronounced with re- 
spect to the gini parameter G, but instead of increas- 
ing linearly through the sample, Figure IT^ suggests that 
Awcm-abs is lowest for galaxies with G < 0.2, peaks for 
intermediate values of G ~ 0.3, and possibly declines 
slightly at the most nucleated values of G > 0.5. This 
high Awcm-abs Sample of high / and moderate - high G 
corresponds reasonably well to a particular morphologi- 
cal sample- the Type 4 galaxies (as identified in §3) which 
occupy a large angular area and tend to have bright neb- 
ulous emission paired with one or more distinct nucle- 
ations. 

These data suggest that larger, more UV luminous (/ 
correlates well with UV luminosity- see §6) galaxies on 
average may have stronger outflows than the rest of the 
population, as might be expected if these galaxies have 
particularly energetic input to their interstellar media 
and are therefore capable of blowing the most energetic 
outflows. We caution however that the major discrepant 
point in both of these trends is for the most nebulous 
sample (i.e. that with both small / and small G) whose 
UV spectra are fainter and typically of slightly lower 
quality. 

Pairing this with the most significant understandable 

We note that the average velocity offset in individual z ~ 3 
LBG spectra with both emission and absorption components is 
typically ~ 650 km s~^- see Shapley et al. (2003) for details. 



20 



Rest-UV Galaxy Morphology 




6500 6600 6500 6600 6500 6600 

Rest Wavelength (A) 

Fig. 16. — Flux-calibrated, composite Ha spectra are plotted 
for each of three bins in the morphological parameters / (size), 
G (gini), 't (multiplicity), and ^ (color dispersion). Spectra are 
numbered according to increasing values of the morphological pa- 
rameter (e.g. G3 is composed of more nucleated sources than Gl). 
The fiducial wavelength of H« and [N ll] A6583Aare indicated by 
dotted lines. Note that the strength of Ho and [Nil] is not signifi- 
cantly different in any of these panels. 
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Fig. 17. — H« fluxes measured for the composite stacks in each of 
the morphological parameters / (size), G (gini), (multiplicity), 
and 5 (color dispersion). Horizontal errorbars represent uncertainty 
in the mean value of a morphological quantity within each bin, 
vertical errorbars represent la uncertainties on the flux based on 
random-draw tests. The dotted line in each panel represents the 
mean Ha flux of the entire sample. 



trend discovered in §5.1 and 5.2, namely the positive 
correlation between G and Lya emission strength, we 
might be led to conclude that Lya and Awcm-abs are 
positively correlated, in constrast to the results of Shap- 
ley et al. (2003) who found that outflow velocities are 
slightly weaker for z ~ 3 LBGs with strong Lya emis- 
sion. This discrepancy illustrates the important point 
that, given the large scatter in all of our correlations, 
correlation is not commutative. That is, while there is 
a general trend that Lya emission strength increases for 
more nucleated objects, this nucleated population is not 
the same as that moderately nucleated, large / popula- 
tion for which At;cm-abs peaks. 

5.4. Rest-Optical Spectroscopic Features 

As part of an ongoing near-IR spectroscopic survey 
(Erb ct al. 2006a) we have obtained rest-frame optical 
spectra in the wavelength regime of Ha and [N ii] for 
19 of the 150 galaxies in the z ~ 2 galaxy sample, and 
use these spectra to explore whether there is any ap- 
parent relation between morphology and Ha flux and/or 
the oxygen abundance as measured by the [Nil] /Ha ra- 
tio (Pettini & Pagel 2004). Given the extremely small 
sample of galaxies with near-IR spectra, we divide these 
19 galaxies into only three bins according to our mor- 
phological parameters and analyze the resulting stacked 
spectra (shown in Figure [TBI) with a method similar to 
that adopted for the rest-UV spectra. 

As indicated by the general similarity of all of the com- 
posites shown in Figure lTBI and plotted more precisely in 
Figure nn we find no significant variation in the strength 
of Ha emission with rest-UV morphology. [N ii] is only 
marginally detected in many of the composite spectra, 
and all variations are well within the uncertainty ex- 
pected based on the noise of the composite spectra. We 
conclude therefore that to within the accuracy permit- 
ted by our small (and hence not fully representative of 



the large distribution of UV morphologies) spectroscopic 
sample the UV morphology of z ^ 2 galaxies is uncorre- 
lated with rest-frame optical spectroscopic features and 
the degree of metal enrichment as parametrized by the 
[Nil] /Ha oxygen abundance estimate. 

6. THE ASSOCIATION OF REST-UV MORPHOLOGIES 
WITH PHOTOMETRICALLY DERIVED PROPERTIES 

Using ground-based UnGTZ and Spitzer- MIPS pho- 
tometry it is possible to calculate the rest-frame, k- 
corrccted luminosities of each of our target galaxies in 
the UV (iieoo) and mid-IR {L^sii) , in addition to es- 
timating the total bolometric luminosity (Lbol) and the 
ratio of IR/UV luminosities (ipiR/iieoo)- The resulting 
spectral energy distribution from these and additional 
JK and Spitzer- IRAC data may then be fit with stellar 
population models to determine the best-fit stellar mass 
(M*), age, and optical extinction {E{B — V)) for a given 
galaxy. Although the risk of confusion is greater in stellar 
population models based on such seeing-limited imaging, 
in almost all cases (except of those of the few most widely 
separated clumps) we find that the isophotes of the tar- 
get galaxies reliably trace the ACS morphology, and all 
components are blurred together into a single object with 
minimal contamination from nearby sources. The com- 
prehensive results of such efforts have been summarized 
by Reddy et al. (2006a); in the present contribution we 
test the degree of association between morphology and 
such photometrically-derived parameters. In this section 
we consider only the 62% of the z ~ 2 sample of galaxies 
for which Spitzer-MIPS 24/im detections directly mea- 
sure the strength of rest-frame mid-IR emission (see fur- 
ther discussion in Reddy et al. 2006b). 

-^^1600 typically ranges from ~ 10^'^ Lq — 1.2 X 10^^ Lq, Ls-g/j 
from ~ 5 X lO^L© - 1.3 X IO^Lq; see Reddy et al. 2006a for 
further details 



Law et al. 2006 



21 



TABLE 2 

Standard deviations from the null hypothesis for 
independence between morphological and photometric 
parameters.^ 



Quantity 


I 


G 




5 


^1600*' 


+6.39 


+4.57 


-0.32 


+3.10 




+0.92 


-0.10 


+0.68 


-0.30 




+0.21 


+0.59 


+0.02 


+0.68 




-3.33 


-1.85 


-0.22 


-1.22 




+0.71 


+0.90 


+0.03 


+0.73 




+0.98 


-0.05 


-0.99 


+0.07 


Age'' 


-0.83 


-0.02 


-0.16 


+1.12 


E{B-VY 


-3.29 


-2.35 


-0.68 


-1.84 



^'+' denotes positive correlation, '-' denotes negative. Bold typeface 
indicate correlations of greater than 90% confidence. 
^Rest-frame ultraviolet luminosity. 
'^Rest-frame mid-IR luminosity (all sources). 
"^Rest-frame mid-IR luminosity (MIPS sources only). 
^Ratio of far-IR to ultraviolet luminosities. 
^Bolomctric luminosity. 
^Stellar mass. 
^Agc of stellar population. 

'Estimated reddening based on far-UV spectral slope. 



We quantify the degree of association using the Spear- 
man non-parametric rank correlation coefhcient r^, val- 
ues of which are assessed in terms of their two-sided sig- 
nificance t. This significance gauges the probability that 
the null hypothesis (i.e. that there is no correlation be- 
tween the two parameters) is true and that any apparent 
correlation is due to random chance alone. In Table |21 we 
calculate t„ for a grid of photometric versus morpho- 
logical parameters (adopting the convention that positive 

indicates a positive correlation, and negative to- a neg- 
ative), highlighting in bold typeface those for which the 
null hypothesis is less than 10% probable (i.e. \t\ < 0.1, 
\t„\ > 1.65). We note from Table H that the null hy- 
pothesis is rejected in very few cases, indicating that in 
general morphology is largely decoupled from photomet- 
rically derived parameters. 

The most significant correlations relate galaxy size 
(/), nucleation (G), and rest-frame ultraviolet luminos- 
ity (ii6oo) in the sense that larger and more strongly 
nucleated galaxies tend to have brighter UV continua. 
In general, these correlations are unsurprising since / ef- 
fectively measures the number of UV bright pixels in a 
galaxy (and a galaxy with large I is hence likely to have 
a large total UV luminosity), and high values of G are 
produced by relatively high concentrations of flux in a 
small percentage of the total number of UV bright pixels 
(which occurs more commonly in galaxies with a high to- 
tal UV luminosity to distribute). However, we note that 
the positive correlation between G and Lieoo may also (at 
least in part) have a more physical explanation; if indeed 
nebulous low G sources are dustier, then UV radiation 
from these galaxies would be more strongly attenuated 
than from less dusty high G galaxies, contributing to the 
observed trend. 

Since t is an awkward statistic to print for extremely small 
values, we instead give results in terms of to- corresponding to the 
number of sigma that t lies out along the wings of a gaussian prob- 
ability distribution. 



The next most significant correlations relate /, G, 
E{B — V) (the estimated reddening based on stellar pop- 
ulation models and the far-UV spectral slope) and the ra- 
tio of IR to UV luminosity (iFiR,/-^i600: a proxy for UV 
attenuation). As quoted in Table |21 / and G are both 
negatively correlated with E{B — V) and iFiR/.^i600 (at 
levels of significance ranging from about 2 — Su) in the 
sense that the most attenuated galaxies tend to have 
smaller UV luminous areas (low I) with weaker nucle- 
ation in the flux distribution (low G) . Although this cor- 
relation contains a great deal of scatter and does not hold 
precisely for every galaxy (i.e. not all "small and neb- 
ulous" galaxies are highly extinguished), it nonetheless 
appears that statistically the most extinguished galax- 
ies tend to be smaller and more nebulous than the gen- 
eral population. As previously hypothesized on the basis 
of Lya emission strength, it may therefore be the case 
that galaxy morphology is dictated in part by the pres- 
ence and location of dust which partially extinguishes the 
UV continuum radiation and causes galaxies to appear 
slightly smaller (to a flxed limiting surface brightness) 
and more nebulous. 

The color dispersion ^ appears to follow the trends 
found for / and G with UV luminosity and extinction (al- 
though at lower significance), possibly due to its known 
positive correlation with G. The multiplicity parame- 
ter 4' however (also loosely correlated with G) exhibits 
no statistically significant associations with any of the 
photometrically derived quantities presented in Table 
121 In addition, we note that contrary to possible ex- 
pectations (e.g that galaxies with higher stellar mass 
might, analagously to local elliptical galaxies, appear 
morphologically more evolved and distinct from those 
with low stellar mass) there is no evidence for a rela- 
tionship between stellar mass and UV morphology, im- 
plying that through a combination of physical processes 
galaxies ranging over two decades in stellar mass from 
Mi, ~ 10^ — 10^"'^ M0 somehow appear visually indistin- 
guishable. However, while no significant overall correla- 
tion was found by the Spearman test, almost all (20/22) 
of the most strongly nucleated sources (G > 0.4) have 
masses greater than 10^*^ M© in contrast to those with 
G < 0.4 which span a full two decades in stellar mass. 

In Table 01 we compare rest-UV morphologies to star 
formation rate (SFR) estimates based on SED fitting, 
UV and IR continuum emission. Although these es- 
timates tend to be loosely correlated with each other, 
there is considerable scatter in individual estimates of 
the SFR for a particular source using different methods. 
The least reliable of our estimates are the SFR estimated 
from normalization of stellar population models (Shap- 
ley et al. 2005; Erb et al. 2006a; Reddy et al. 2006a), 
for which we find no significant correlations with mor- 
phological parameters (due either to the intrinsic lack of 
such a correlation, or to the uncertainty inherent in the 
SFR estimate given by this method). The next most re- 
liable estimate is that based on UV continuum flux, for 
which we calculate the uncorrected star formation rate 
based on UnGTZ photometry and the conversion factor 
given by Kennicutt (1998). As shown in Tabled the 
uncorrected UV SFR is strongly correlated with all mor- 
phological parameters save the multiplicity ^E*. However, 
this is not surprising given the strong association between 
total UV luminosity and morphology shown in Table 12 
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TABLE 3 

Standard deviations from the null hypothesis for 
independence of morphological parameters and star 
formation.^ 



Quantity 


/ 


G 


-I- 




SFRsed'' 


-0.32 


+0.59 


-0.15 


+0.30 


SFRuv" 


+8.04 


+4.77 


+0.46 


+3.15 


SFRiRd 


+0.21 


+0.59 


+0.02 


+0.47 


SFRuv+ir'' 


+0.92 


+1.05 


+0.03 


+0.87 


SFRuv+iR^ 


+5.67 


+2.17 


+0.80 


+ 1.76 




+1.24 


+1.62 


+0.88 


+1.01 




-1.12 


-0.43 


-0.04 


-0.70 



^ + ' denotes positive correlation, '-' denotes negative. 

^Star formation rate estimated from SED fitting. 

'^Star formation rate estimated from UV photometry. 

^Star formation rate estimated from IR photometry. 

^Star formation rate estimated from UV+IR photometry, MIPS de- 
tected sources only. 

^Star formation rate estimated from UV+IR photometry, MIPS un- 
detected sources only. 

^Specific star formation rate, cj} — SFRuv+ir/A^,^ . MIPS detected 
sources only. 

^^Spccific star formation rate, (f> — SFRuv+ir/A/j, . MIPS undetected 
sources only. 



While it is likely that large (high /) , nucleated (high G) 
objects have a higher than average rate of uncorrected 
star formation it is not necessarily true that the same 
holds for the total star formation rate. Indeed, in the 
majority of cases obscured star formation comprises the 
bulk of the total. Rather than assuming an average ex- 
tinction factor to estimate this from the UV luminosity, 
wc estimate the obscured SFR from the IR luminosity 
LiR calculated by Reddy et al. (2006b), adopting the 
LiR— SFR conversion factor given by Kennicutt (1998). 
As shown in TableO obscured star formation has far less 
compelling associations with UV morphology- none are 
found to have statistical significance. 

The most accurate estimate of the total star formation 
rate may be found by summing the rates determined 
from UV and IR luminosities. We break this compar- 
ison into two cases- galaxies for which the IR-derived 
SFR has been directly measured using Spitzer-MIPS de- 
tections at 24/im, and galaxies which are undetected at 
24/im and an estimate of the IR-derived SFR has been 
computed based upon the average bolometric luminos- 
ity of MlPS-undetected sources" (Lbol = 6 x lO^'^L©, 
Reddy et al. 2006b). For the first of these cases, no 
significant correlations are found with either the total 
or specific (stellar mass-normalized) star formation rate. 
We note however that the specific star formation rate is 
on the threshhold of statistical significance for a positive 
correlation with G. If the eleven galaxies with stellar 
populations younger than 40 Myr (i.e. for which stel- 
lar population fits are the most uncertain) are neglected 
in our analysis, this correlation becomes statistically sig- 
nificant, deviating from the null hypothesis by approxi- 
mately 2.1(7. While therefore there is mild evidence for a 
positive correlation between UV nucleation and net star 
formation for galaxies directly detected at 24/im, this 
association is not at present statistically compelling. In 

LjR effectively scales with Lbol- 



contrast, for the 24^m- undetected sources all of the asso- 
ciations between morphology and uncorrected UV SFR 
are recovered for the total estimated SFR (although this 
trend is smeared out by the division by mass for the 
specific SFR), largely because the obscured SFR extrap- 
olated from total bolometric luminosity contributes only 
marginally to the total in these cases. 

7. A COMPARISON OF MORPHOLOGIES WITH OTHER 
GALAXY SAMPLES 

7.L AGN/QSO 

As given in the table of AGN/QSO in GOODS-N pre- 
sented by Reddy et al. (2006a, their Table 3), there are 
a total of nine spectroscopically confirmed AGN/QSO in 
our sample in the redshift range z — 1.8—3.4, six of which 
are directly detected in X-ray emission (although one 
very weakly) and three of which have no X-ray counter- 
part to a depth of 2 Ms (based on the catalog of Alexan- 
der et al. 2003) but are confirmed AGN based upon 
high-ionization optical emission lines and/or power-law 
mid-IR SEDs. As noted by Reddy et al. (2006a), the 
morphologies of the X-ray undetected sample are more 
disturbed than their directly detected counterparts. As 
we show in Figure ^1 all three X-ray undetected sources 
have G < 0.4, while the the five detected sources all have 
strongly nucleated values G > 0.38 (four with G > 0.47) 
and the two QSO in the sample (BMZ 1083 and MD 39) 
have the highest values of any object considered in this 
contribution at G = 0.63 and 0.88 respectively (close to 
the stellarity limit of G ^ 0.90). This suggests a possible 
correlation between X-ray luminosity and UV nucleation, 
in the sense that AGN which produce obvious X-ray radi- 
ation also contribute sufficiently to the total UV output 
of their host galaxies that the centralized AGN radia- 
tion visibly affects the apparent nucleation of the UV 
light profile. Given the extremely small size of our sam- 
ple, it is difficult to asses the global applicability of this 
correlation although we note for comparison that in a 
larger sample of 31 AGN in the redshift range z ^ 2 — A 
Akiyama (2005; their Fig. 3) found that AGN with the 
most centrally concentrated light distributions tended to 
be X-ray bright (L2-iokeV > 10'*'' erg s~*) broad-hne 
sources, while the population of X-ray bright narrow- 
line sources and X-ray faint sources tended to be slightly 
less centrally concentrated. Likewise, recent studies of 
AGN at lower redshifts (z < 1.3) by Pierce et al. (2006) 
and Grogin et al. (2005) also found the rest-frame op- 
tical morphologies of X-ray luminous AGN to be more 
compact than those of IR-selected AGN or ordinary field 
galaxies. 

Not only are the X-ray undetected AGN more nebu- 
lous than their X-ray bright counterparts, we find that 
they also appear to have a greater number of compo- 
nents to their UV light, exhibiting a range of multiplic- 
ities ^I* — 4.6 — 11.8 as compared to the ^' < 1.8 found 
for all directly detected X-ray sources. These irregular 
morphologies, coupled with a power-law SED longward 
of 3 /im (indicating the presence of warm dust) yet a 
UV SED well-fit by a simple stellar population suggest 
(Reddy et al. 2006b) that these sources may be obscured 
AGN whose UV emission is dominated by spatially ex- 
tended star formation rather than a central active nu- 
cleus. 
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Fig. 18. — Histrogram showing the relative UV nucleation 
for the z = 1.8 — 3.4 galaxy sample (red histogram) versus X- 
ray detected (green histogram) and X-ray undetected AGN/QSO 
(blue histogram). The nucleation of the X-ray undetected sam- 
ple is broadly consistent with the median nucleation of the overall 
z = 1.8 — 3.4 population, while X-ray detected galaxies are more 
strongly nucleated. 



7.2. IR-selected BzK Galaxies 

One well-studied near-IR color selected galaxy sam- 
ple is the BzK- selected catalog (Daddi et al. 2004), 
for which many sources (fifty two sources brighter than 
K — 2\ \ii the redshift range z = 1.8 — 2.6) simulta- 
neously satisfy both the BzK selection criteria and our 
optical UnGTZ color selection criteria. In Figure ITOI we 
show the morphologies of those GOODS-N BzK galaxies 
which are comparably bright in K band and and have 
well-determined photometric redshifts, yet fail to meet 
the optical selection criteria. Qualitatively, we note that 
the BzK galaxies not in the optically-selected sample 
appear morphologically very similar to those in Figure 
12 as might be expected given the large general overlap 
between the two samples and suggesting that these few 
galaxies may fall just outside of the optical selection cri- 
teria. 

We quantify these morphological differences in Figure 
1201 using photometric redshift estimates to scale the sur- 
face brightness selection algorithm. As expected based 
on Figure 1191 the average morphology of BzK sources 
which are also UnGTZ sources (filled triangles) is identi- 
cal to that of the overall UnGTZ sample (filled squares) 
to within the uncertainties in the mean. The sample of 
eleven BzK sources which do not meet the UnGTZ criteria 
(empty triangles) has slightly fewer high-surface bright- 
ness pixels (lower /) and is slightly less nucleated (lower 
G) than the other two samples, but at a confidence level 
of only 1 — 2a. The apparent multiplicity 5" and color 
dispersion ^ remain approximately constant among all of 
these samples to within the uncertainty. 

7.3. IR-selected Distant Red Galaxies 

Distant red galaxies (DRG's, Franx et al. 2003) con- 
stitute another major IR-selected sample of high red- 
shift galaxies for which there are five sources brighter 



than K = 21 'm the redshift range z — 1.8 — 2.6 which 
meet the UnGTZ selection criteria and have secure spec- 
troscopic redshifts, and eleven sources which do not but 
for which we measure reliable photometric redshifts. As 
illustrated by Figure [21 the DRG population is fainter 
and more nebulous in UV emission than any of the other 
galaxy samples considered, consistent with what may be 
expected for particularly dusty IR-bright galaxies given 
the apparent correlation between UV nucleation and dust 
extinction. 

As for the BzK sample, the overlapping sample of 
C/„G7?.-selected DRGs has morphological coefficients sim- 
ilar to the bulk of the UnGTZ population to within l — 2a 
(although it is on average slightly slightly smaller and 
more nebulous, see Fig. I20|l . It is those DRG not selected 
by UnGTZ criteria however which particularly stand out- 
these galaxies (Fig. (201 open circles) have much smaller 
/ and G to a high level of confidence and in some cases 
(e.g. DRG 14) are barely detected in our HST-ACS imag- 
ing data. Largely, this difference is a reflection of the 
much fainter UV luminosity of the DRG sample- typi- 
cal DRG which do not meet the optical selection criteria 
have UV luminosities Liqqq < (and commonly 

^1600 < 10^X0), as compared to the median luminosity 
for the optically selected sample iieoo ~ 5 x 10^'-' Lq. 
Similarly to the BzK sample, the apparent multiplicity 

for DRGs is again consistent with the UnGTZ sample, 
although the color dispersion ^ is roughly 2<7 lower and 
is consistent with zero (the value expected for a pure 
measurement of the background sky). 

7.4. Submillimeter Galaxies 

The submillimeter galaxy (SMG) population offers the 
opportunity to compare the rest-UV morphologies of 
optically-selected galaxies with those in a similar redshift 
range (z ^ 2 — 3) which are selected on the basis of sub- 
millimeter flux. Using the catalog of SMG coordinates 
and spectroscopic redshifts compiled by Chapman et al. 
(2005), we apply our morphological analysis to a sam- 
ple of these galaxies in the GOODS-N field. The wide 
range of morphological types covered by the SMG selec- 
tion criteria are shown in Figure l^ and range from single, 
nucleated sources to extremely faint and nebulous, and 
include one strong high-redshift spiral galaxy candidate 
(BZK 294- see Dawson et al. 2003). On average however, 
the SMG population has a morphology similar to that of 
the [/„G7?.-selected sample (Fig. I20|) with values of 4' and 
^ consistent to within the uncertainty and apparent sizes 
/ and nucleations G only slightly smaller and more neb- 
ulous than the [/„G7^-selected sample. These results are 
consistent (since our size parameter / is closely related 
to surface brightness) with the findings of Chapman et 
al. (2003) and Conselice et al. (2003b), whose 775T-STIS 
imaging led them to conclude that the SMG sample had 
lower surface brightness than typical z ~ 3 LBGs, yet 
was slightly larger for the reduced surface brightness. 

Additionally, Conselice et al. (2003b) also find that 
the asymmetry index A of their sample of eleven SMGs is 
slightly greater on average than that of optically-selected 
z ^ ?> LBGs, which these authors conclude implies 
a greater major merger fraction for the submillimeter- 
bright sample. However, we find no statistically sig- 
nificant change in our multipliticy parameter 5" (which 
most closely measures the morphological irregularity of a 
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Fig. 19. — As Figure l^l but for IR-sclcctcd BzK sources brighter than A' = 21 which arc not selected by the rest-UV color selection 
criteria. Redshifts given are photometric. 
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Fig. 20. — Comparative average size (/), nucleation (G), multiplicity (']/), and color dispersion (g) for optically selected UnGTZ, BzK, 
DRG, and SMG galaxy samples. Filled squares represent all z 2 GT^-selected galaxies, filled triangles and circles the respective subsets 
which also fulfill BzK and DRG criteria, open triangles and circles galaxies which fulfill the BzK and DRG criteria (respectively) which do 
not meet the UnGTZ selection criteria, and open squares all galaxies which fulfill the SMG critera. All galaxies are limited to the redshift 
range z = 1.8 — 2.6 (using photometric redshifts where spectroscopic are unavailable), and the BzK and DRG galaxy samples are limited 
those with no x-ray counterparts to 2 Ms and K magnitudes K < 21. Errorbars indicate the la uncertainty in the mean values for each 
population. 



source) . If the morphologies of the SMG sample are gov- 
erned by the same physical processes as the optically- 
selected galaxies (as seems likely given the two sam- 
ples appear to fall on the same trend of dustiness ver- 
sus bolometric luminosity for star-forming galaxies at 
z ~ 2; Reddy et al. 2006b) then when considered non- 
parametrically we conclude that the SMG sample may be 
no more likely to be dominated by major mergers than 
other z 2 — 3 galaxy samples. 

8. DISCUSSION 

Through a comprehensive statistical analysis, we have 
found two key trends with morphology. First, the de- 
gree of nucleation/nebulosity is related to the UV/optical 
extinction of the source (as parametrized through both 
the value of E{B — V) calculated from stellar population 
models and the ratio of IR to UV luminosity) , suggest- 
ing that more nebulous sources may appear as such since 
their greater quantities of dust obscure a great quantity 
of the UV light. Second, physically larger, more UV 



luminous sources have rest-UV emission and absorption 
lines seperated by a greater velocity, indicating that they 
may drive more energetic outflows than their lower lumi- 
nosity counterparts. However, we note that these trends 
are typically on the order of 3 — 4fT significance, and 
are far from being unambiguously convincing morpholog- 
ical differences. Indeed, given the relatively large sam- 
ple of galaxies and volume of spectrophotometric data 
compiled to date, the overall lack of correlation between 
morphology and fundamental galactic properties such as 
stellar mass, SFR, and outflows is more remarkable. To 
some degree this may be due to the unknown distribu- 
tion of viewing angles, but whether as a result of this 
complication or a more fundamental process the rest- 
UV morphologies of high redshift galaxies generally do 
not appear to contain a great deal of separable informa- 
tion. If morphologies truly mean so little, this may help 
explain the lack of correlation seen between elongated 
morphologies and major axis velocity shear by Erb et al. 
(2004), who found that morphologically elongated galax- 
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Fig. 21. — As Figurc|^ but for IR-selected DRG sources brighter than K = 21. Top row: DRG which are also selected by rest-UV color 
selection criteria, Bottom rows: DRG which are not selected by rest-UV color selection criteria. Redshifts given for galaxies not selected 
by the rest-UV criteria are photometric. 
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Fig. 22.— As Figure^ but for all SMGs on the HST-ACS image frame of the GOODS-N field in the rcdshift range z = 1.8 - 2.6 (plus 
one source at z = 2.91). Positions and spectroscopic redshifts are drawn from the catalog of Chapman et al. (2005). Galaxies which also 
fulfill UnGTZ, BzK, or DRG selection criteria are labelled accordingly. Note that the lower right object in this figure (BM 1083) is a known 
QSO. 



ies are no more likely than compact galaxies to exhibit 
kinematic signatures of rotation. 

It is interesting to note the implications of these find- 
ings for the major-merger hypothesis which posits that, 
similarly to luminous and ultraluminous infrared galax- 
ies (LIRGs and ULIRGs) in the local universe, particu- 
larly irregular morphologies correspond to major galaxy 
mergers (e.g. Conselice et al. 2003, Chapman et al. 
2003, LPM04, Lotz et al. 2006). In particular, Conselice 
et al. (2003a) adopt the asymmetry selection criterion 
A > 0.35 to identify major mergers, and interpret the 
increase in galaxy irregularity with redshift out to z ~ 3 
to indicate a higher merger fraction in the early universe 
and evidence for hierarchical formation. If the most ir- 
regular galaxies truly represent such major events how- 
ever, we might reasonably expect these galaxies to dif- 
fer significantly from their non-merging counterparts in 
some additional way, perhaps by displaying tidally en- 
hanced rates of star formation or broadened interstellar 



absorption lines due to the juxtaposition of the ISM of 
two or more galaxies with relative velocity differences of 
order a few hundred km s^^. However, we find no ev- 
idence for such a distinction between any of a variety 
of morphological samples, suggesting that either major 
mergers at z ~ 2 — 3 are largely indistinguishable from 
the non-merging sample, that aZ/ z ~ 2 — 3 galaxies are 
experiencing major merger induced starbursts, or simply 
that irregular morphology (at least to within ^ 13 kpc at 
rest-UV wavelengths) is not a reliable indicator of a ma- 
jor merger. While it is possible that major mergers may 
be distinguishable on the basis of multiple components 
beyond our 13 kpc detection threshhold, it seems unlikely 
that such systems would produce changes as profound as 
may be found on closer passes. 

Supposing then that irregular, multi-component mor- 
phologies do not represent bursts of star formation in 
the orbiting galaxies and satellite galaxies of a major 
merger, another possibility is that these systems may in- 
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stead represent patchy star formation occuring from the 
collapse of local instabilities within galactic molecular 
clouds. Combined with the complicating effects of dust 
and viewing angle, such an explanation may fit the lo- 
cal LBG-analog VV If 4, based on the UV and near-IR 
imaging of Goldader et al. (2002). ft is also possible 
that, in some small number of cases, multi-component 
morphologies may reflect objects which are nearby in 
projection, but unrelated due to large differences in red- 
shift. While individual components of multi-component 
objects have not been spectroscopically confirmed to be 
associated however, the angular distribution of sources 
suggests that such projection effects should contribute 
minimally to the total population. 

Finally, we note that while the rest-UV and rest-optical 
morphologies of z ~ 2 — 3 galaxies are similar in the 
majority of cases (as shown by Dickinson 2000) and ap- 
pear to contain little separable information, there is a 
small sub-sample of these galaxies for which the opti- 
cal morphologies appear considerably more regular and 



"evolved" than their UV counterparts (Toft et al. 2005). 
It may be the case that this small population of galaxies 
has rest-optical morphologies which are easier to inter- 
pret than the majority of the galaxy sample, and may 
have useful correlations between rest-optical morphology 
and spectrophotometric or kinematic indices. In particu- 
lar, such galaxies are intriguing possible targets for future 
kinematic study with the aid of AO-assisted integral-field 
spectroscopy (e.g. Forster-Schreiber et al. 2006; Law et 
al. 2006), which could demonstrate whether these galax- 
ies are any more likely than the rest of the galaxy sample 
to contain tractable meaning in their luminous spatial 
structures. 
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